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Abstract. We develop a theory of particle acceleration inside relativistic rotating electron-positron force-free 
jets with spiral magnetic fields. We considered perturbation of the stationary magnetic field structure and found 
| that acceleration takes place in the regions where the Alfven resonant condition with the eigenmodes in the jet is 

t-H ■ fulfilled, i.e. where the local Alfven speed is equal to the phase speed of an eigenmode. The acceleration mechanism 

is regular acceleration by the electric field of the electromagnetic wave, which is the eigenmode of the force-free 
cylindrical jet: particles drift out of the region of the large wave amplitude near the Alfven resonant surface and 
gain energy. Acceleration in the strong electric field near the Alfven resonance and synchrotron losses combine 
to form a power-law energy spectrum of ultra-relativistic electrons and positrons with index between 2 and 3 
depending upon the initial energy of the injected particles. The power law distribution ranges from ~ lOMeV to 
Q_,' ~ 10 3 MeV. 
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1. Introduction 

The nature of the Active Galactic Nuclei (AGN) is still unclear. The most common point of view is that there is a 
supermassive black hole in the centre of the active galaxy with a mass approximately 10 s to 10 9 Mq (Rees, 1984). 
Accretion onto the black hole drives powerful turbulent processes in accretion discs which effectively heat the plasma 
and generate a magnetic field of the order of 10 4 G. Ejection of plasma from discs with frozen magnetic field and 
radiation from the disc lead to the formation of collimated streams which transfer energy to large distances, up to 
10 5 pc. In addition, the processes of generation of electron-positron plasma may occur in the magnetospheres of black 
holes as well as in the vicinity of accretion discs. Particles are created in collisions of high-energy photons produced 
by inverse Compton scattering of ultraviolet photons emitted from the disc with fast particles accelerated by magnetic 
reconnection within the disc or within the magnetosphere of the rotating black hole, where the equilibrium charge 
density changes its sign (Beskin et al. 1992; Hirotani & Okamoto, 1998). Electron-positron plasma forms pinched 
streams (jets) which have radii of the order of one parsec. The electromagnetic energy density in these jets greatly 
exceeds the total energy density of electrons and positrons, including the rest energy. The angular momentum from 
the black hole and the accretion disc is transmitted to the jet via magnetic coupling, where a strong radial electric 
field E r arises. Drift in this field and in the longitudinal magnetic field, B z , frozen into plasma leads to differential 
rotation. The jet transfers an electric current producing an azimuthal magnetic field, B^, and magnetic field lines twist 
into spirals. The energy in the jet is transferred by the Pointing vector which is proportional to the product of the 
radial electric field and the azimuthal magnetic field. Therefore, there is no problem of energy transfer along the jet 
to large distances compared to the case where the main energy is in the hydrodynamic motion. 

The amazing facts about the jets from AGN are their high collimation (the ratio of the length to the diameter is a 
few orders of magnitude), the stability, and apparent superluminal velocities of distinct knots (Begelman et al. 1984). 
The hydrodynamic and MHD stability of jets was investigated in many studies (Blandford & Pringle, 1976; Torricclli- 
Ciamponi & Petrini, 1990; Appl & Camenzind, 1992). The stability of the electron-positron jets, rotating and moving 
with relativistic speed and surrounded by a dense medium was studied by Istomin & Pariev (1994) and Istomin & 
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Pariev (1996). Such jets were shown to be stable with respect to axisymmetric as well as spiral perturbations. The 
physical reason for the stability is the shear of the magnetic field, which in the case of low hydrodynamic pressure 
stabilizes small oscillations. The density of the media surrounding the jet, p, should satisfy the condition p >• B 2 /(Aire 2 ) 
for the inertia of surrounding gas to stabilize the relativistic jet. On parsec scales, a typical value of B ~ 10~ 2 G gives 
a proton density n >• 0.05 cm~ 3 . The characteristic density of gas in galactic nuclei is n ~ 10 cm -3 . Therefore, the 
approximation of stationary jet walls (Istomin & Pariev, 1994; 1996) is well justified. A relativistic force-free jet with 
infinite radial extent or with vacuum outside can develop weak instabilities (Lyubarskii, 1999; Li, 2000). In a stable 
jet, perturbations do not increase with time (Imw = 0) or have a small decay rate (Imw w 10~ 2 Rew) because of 
the resonance with Alfven waves u)' = ftf, c. When the specific energy density and pressure of the plasma are much 
less than the energy density of the magnetic field, the Alfven velocity is equal to the speed of light in vacuum c [ui 1 
and k' are the frequency and the wave vector in the plasma rest frame). The resonance condition is fulfilled on the 
specific magnetic surface and, in the case of cylindrical jet, at a definite distance from the axis. In the vicinity of that 
surface the magnetic and electric fields of the wave are large. Particles in the jet are accelerated by the electric field 
there, absorbing the energy of the perturbation. Thus, the stability of the jet is directly related to the production of 
energetic particles in the jet. This resolves the problem of in situ acceleration of electrons and positrons producing 
synchrotron emission from knots observed along jets. In fact, as is known (see Begelman at al., 1984), the energetic 
particles accelerated in the central region cannot survive far along the jet; they must lose their transverse momentum by 
synchrotron radiation already at the basement of the jet. The time for an electron to lose its perpendicular momentum 
due to synchrotron radiation is given by 

_ 3 gyy 

s — A 2 

4 r c uj z m 

Here 7|| is the longitudinal Lorentz-factor of the particle, 7|| = (1 — f|/c 2 ) -1 / 2 ; r c is the classical radius of the electron 
r c = e 2 /m c c 2 ; ujbo = eB/m c c is the electron cyclotron frequency. Given 7|| = 10, B — 10 2 G, r s is approximately 
3 • 10 5 sec, whereas the typical time of the flyby through the region of strong fields T{ s=s £/c rs 10 6 sec. Here I 
is the diameter of the jet, t s=s 3 • 10 16 cm. Qualitatively the possibility of particle acceleration by electromagnetic 
disturbances propagating along a jet was proposed by Bisnovatyi-Kogan & Lovelace (1995). However, the authors of 
this work needed an ad hoc vacuum gap between conducting jet and surrounding medium, where the acceleration is 
possible due to the presence of the component of electric field parallel to the magnetic field. 

In this paper we assume that the energy transferred along the jet is mainly the electromagnetic energy which can 
be transmitted for long distances along the jets as if along wires. In this case it is very possible that part of that 
energy is in the waves propagating along the jet, which are eigenmodes of a cylindrical beam. The source of the wave 
motion consists of non-stationary processes in the magnetospheres of the black hole and the accretion disc. Short time 
variability on scales from days to months is actually observed in AGN (Mushotzky et al. 1993; Witzel et al., 1993). 
Such a picture also gives a natural explanation of the supcrluminal motions but requires neither special orientation of 
the jet nor a relativistic jet with high speed (/3 = u/c > \/2/2). As was shown in papers by Istomin & Pariev (1994) 
and Istomin & Pariev (1996) there exist so called standing eigenmodes in jets (v gr0 up = 0) which do not propagate 
with finite velocity along the jet but are only subject to diffuse spreading and, therefore, their amplitude is maximal. 
The phase velocity of these modes is greater than c. Wave crests move along the jet with superluminal velocity causing 
the acceleration of the particles on the Alfven surface. These regions might be the observed bright knots with typical 
sizes of the order of the wavelength of the standing wave, which is about the radius of the jet. 

There is some observational evidence in the case of 3C273, M87, and 3C345 that their jets are electron-positron 
rather than electron-proton (Morrison & Sadun, 1992; Reynolds et al., 1996; Hirotani et al., 2000). The recent ob- 
servation of circular polarization of radio emission from 3C279, 3C273, 3C84 and PKS0528+134 (Wardle et al., 1998 
) supports the same conclusion. However, there are also indications that jets in Optically Violently Variable quasars 
cannot consist solely of e + e _ pairs because they would produce much larger soft X-ray luminosities than are observed. 
On the other hand, models with jets consisting solely of proton-electron plasma are excluded, since they predict much 
weaker nonthermal X-radiation than observed in Optically Violently Variable quasars (Sikora & Madejski, 2000). 

We consider the acceleration processes near the Alfven resonant surface and calculate the spectrum of accelerated 
particles. In an accompanying paper (Pariev et al. 2002), we calculate the polarization of synchrotron radiation 
generated by these particles. 

2. Acceleration of particles 

The equations of motion near the Alfven surface are rather complicated. Nevertheless, we can use drift equa- 
tions (Sivukhin, 1965) because the Larmor radius of relativistic electrons and positrons r c = £/eB w 10 6 cm (if 
we take B w 10~ 2 and 7 w 10) is much smaller than radius of the jet and the width of the resonant surface ro (see 
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expression (28)): 

| - if + £!■ >< B] + ||[B x (BV)B] + J^L[B x VB]; 
d£\ _ dr £v\ dB 

Here the electric field, E, is assumed to be small compared with the magnetic field, B, and £ is the energy of the 
particle. The index e.m. means the energy change due to the action of the electromagnetic fields, E and B on the 
particle. We did not write the third equation here, the conservation of the adiabatic invariant J± = p\/B, which 
is usually used as one of the drift equations. The adiabatic invariant is not conserved because of synchrotron losses. 
But the characteristic time of the change of the adiabatic invariant is much larger than the period of the cyclotron 
rotation and we can use the drift approximation of the motion of the particle as a result of the separation of the motion 
into fast cyclotron rotation and slow motion of the guiding centre with subsequent averaging over the fast rotation. 
The conditions for using the drift approximation are a small change of the cyclotron frequency ujb during cyclotron 
period 2tt/ujb and also a small change of the particle cyclotron radius on the scale of the cyclotron radius (Alfven & 
Falthammar, 1963): 

2tt duiB I ___ | ,„ . 

^T<^b; r B \Vr B \<^r B . (la) 

lob at 

A rigorous derivation of equations (1) under conditions (la) was given by Sivukhin (1965) and Morozov & Solov'ev 
(1966). Let us demonstrate that the second equation (1) does not require the conservation of the adiabatic invariant 
J±. The particle velocity is split into two parts v = dr/dt + v c . Here dr/dt is the velocity of the guiding centre and 
v c is the velocity of cyclotron rotation, v c = v±. 

( d£\ „ dr „ dr euos f^„ dr „ eoJs f , 

— = evE = e— E + ev c E = e— E H & Edl = e— E H / curl Eds. 

V dt J c m dt dt 2tt J dt 2tt J 

Using the Maxwell equation, curlE = — (<9B/dt)/c, and the relation s = — tt(v^/ujb) 2 B/B we obtain the second 
equation of the drift approximation (1). The first equation of system (1) is the equation of the motion of the guiding 
centre dr/dt and contains the motion along the magnetic field line, the electric drift, gradient and centrifugal drifts. 
All drifts and the equation of energy change are independent of the conservation of the adiabatic invariant J±. The 
value of J± is conserved in the slowly varying electromagnetic fields B(r, t), E(r,i), but is not conserved under the 
action of additional forces including radiative forces. The small change of the adiabatic invariant during the period of 
the cyclotron rotation does not violate the condition of applicability of the drift approximation (la). 

The electromagnetic fields are equal to the sum of the stationary fields (without subscript) and the wave fields 
(subscript 1). For a cylindrical jet the stationary configuration of fields is (Istomin & Pariev, 1994) (c=l): 

B = B z e z + B e ; B<f, = Q F rB z : 

u = KB + fi F re ; (2) 

E = -fi F r[e^xB], 

where r, <j>, z are cylindrical coordinates. In components 

B = (0,fi F r,l)B,, 

E = (-n F r,0,0)B z . 

Here fl F and K are functions of r, and B z does not depend on r. The stationary electric field, which has an absolute 
value of Q F rB z , is not small compared to the magnetic field when Q F r <~ 1; therefore we must consider our equations 
in the frame moving with the plasma where E = 0. There are many such reference frames because there is an arbitrary 
parameter K(r) in the expression for the plasma velocity which determines the radial profile of longitudinal velocity. 
K(r)B z is the velocity along the magnetic field which is not related to the rotation. We choose, however, the velocity 
u which minimizes the kinetic energy of the plasma in the stationary reference frame. 

_ (0,» F r,-(» F r) 2 ) 
U ~ 1 + (fFr) 2 ' [6 > 

It is probable that the plasma moves with that velocity in real jets (Frolov & Novikov 1998). THE fields in the wave 
calculated in terms of the Lagrangian radial displacement £ are (Istomin & Pariev 1996) (here we drop out the phase 
coefficient exp(— iut + ikz + imcf))): 

B rl = iB z F£; 



4 



A.R. Beresnyak et al.: Relativistic parsec-scale jets: I 



B^ = -B z (n F rf+^(n F r) + ^D); 

B rl = B,(S! p r| + fi(^r)-|D): 

i?0i = —iB z {u> — mf2 F )£; 
= iB z QFr{uj + k)£; 

where 

Z ? = r |(0 F ( w + fc )--)- e (n F (a; + fc ) + 5), 

f = k + m n F , 

S = uj 2 -k 2 - m 2 /r 2 . 

For Lagrangian displacement £ and dimensionless pressure disturbance = 47rPi / £? 2 we have a system of differential 
equations obtained by Istomin & Pariev (1996): 

Here 

Ci = ^ (mO F - (f7 F r) 2 (cj + fc)); 

o> 2 - fc 2 - m 2 /r 2 
Go = — - 



+ k 

C 3 = -(uj + k)(A 2 -m F2 ); 
A = k-uj + 2mfl F - (Q F r) 2 (uj + k). 
Let us expand Eqs. (5) near the Alfven surface, i.e. near the point ta where A(rpA = 0. Denoting x = r — r& we get 



. The general solution of these equations is 



/ C 2 \nx \ (C 2 , 



where A' = 

dr 



pj \Ci\nx-A'J *\a 

where a\ and a 2 are constants. They are not arbitrary; rather they are fixed by the boundary conditions on £ and 
obtained from solving general equation (5). We see that £ and p* have a logarithmic singularity at the Alfven point. 
Later we adhere to our expansion over r— r\, assuming r = r\, A = whenever the quantities under consideration have 
no singularities at r = r\. Let us calculate dr/dt and (d£ / dt) e . m . to first order in the fields of the wave. Remembering 
that E = in our reference frame and knowing that terms with the energy £ contain a small parameter, either the 
ratio of the Larmor radius to the jet radius or to the characteristic length of variations of zero order fields, we can 
simplify the drift equations as follows, 

'dr' 



dt 



= v { \e z 

o 



dr\ Bi B(B • Bj) [Ei x Bl 



The subscript denotes the zero-order term and 1 denotes the first order term. The force of inertia is also present in 
the expression of dr/dt (see Sivukhin, 1965). The term with inertia force contains the electron mass and has the same 
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order as terms with energy and is omitted. It is clear that in the force-free approximation the impact of the inertia 
force to the trajectory is negligible. When we substitute Bi and Ei into Eq. (8) we get 

-A =h-D^ (*) 



where 7 = \/\ + (f2 F r) 2 is the Lorentz factor corresponding to the velocity of the moving frame, £ depends on z, </>, t 
by exp{— iujt + ikz + im<f)}. Using the Lorentz transformation, the phase of the wave \t "seen" by the particle is found 

F F F 
to be \f = t' -\ z'. Assuming dz/dt = vn we get * = (v\\ — l)—t'. As we see, the coefficient for the phase coincides 

with the coefficient in Eq. (8*). This confirms our calculations because £ is the Lagrangian displacement. The question 
of the trajectory of the particle will be considered later. 
Now we proceed with calculating (d£ / dt) e . m .: 

d£\ „ fdr\ £v 2 , dB x 

After substitution of expressions for fields into (9) we get 

f) =^i(-fU(. + (10) 



or, combining with (8*) 



d£\ = _£vi_ nu, + k)[F-u>Y ^ (U) 



The inertia force gives the term without the coefficient (u> + k)[F— u>]/S. This coefficient is large if the fast magnetosonic 
resonance surface S = happens to be close to the Alfven resonance r^. Our calculations (Istomin & Pariev 1996) 
show that S is indeed small at r — when the damping rate of the eigenmodes is maximal. The maximal damping 
rate is still small compared to the frequency u. Rapid damping of energy in these modes leads to a high efficiency of 
particle acceleration, which depends, as we will see later, on the smallness of S (Q cx 1/S in Eq. (23) must be large in 
order to get high particle acceleration). The inertia force is of the order of £/r, which may not give large acceleration 
because Ar << r^. As we see, £ is proportional to r if we do not take into account the dependence of the right hand 
side on £. (Note again that we consider all values at the Alfven resonance, r = ta-) 

Therefore it makes sense to consider the motion of the particle averaged over many periods. Taking the r-component 
of (8*), we have 

dr * F 

-= iuJ *Z(x)e™ t , W * = ( V ||-1)-. (12) 

Near the Alfven resonance £(x) has a logarithmic singularity (Eq. (7)); therefore, Eq. (12) can be written as 

dx 

— =w*[Ailnxsin(w*f) + 5iCos(w*f)], (13) 

where u>*A\ and ui*B\ are the dimcnsionless amplitudes of the perturbation (remember that c = 1). Let us average 
Eq. (13), expanding it in small amplitude of oscillations of a particle in the field of the wave. We denote quantities of 
zero, first, second order as Xo, x\, X2- 

= w*L4i \nx sm(oj*t) + B\ cos(a>*t)], 
X\ = —Ai lnx cos(uj*t) + Bi sm(u*t), 

dx2, 

—jj- = uj*[Ai In \xq — AItixq cos(uj*t) + B\ sm{u*t)\ sin(w*t) + B\ cos(w*t)], 



dX2 _ * ^ 1 
dt 1 1 2x ' 



dx 

Now let us replace x$ by x and — - — by — to obtain a smoothed equation of motion 

dt ' dt 

dx _ ^AiBi 
~dl~ U 2x ' 
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which has the solution 

x = y/^AxBtf. (14) 

The result for the smoothed particle motion (14) was checked by the direct numerical integration of Eq. (13). 

From Eqs. (11), (14) it follows that the particle gains energy drifting along the r axis. The rate of acceleration 
is proportional to the transverse momentum. Because of synchrotron losses, p± and J± are monotonically decreasing 
according to (u>b = ^Bo m cC 2 /£) 



dJ ± 



1 - 



J±, 



(15) 



this leads to pu > p±. Such an anisotropic distribution can become unstable to the excitation of electromagnetic waves 
which results in isotropization of the distribution function in momentum space. The synchrotron losses and anisotropic 
instabilities lead to the violation of the adiabatic invariant conservation J± = constant. 

First of all we consider instabilities in a magnetically dominated plasma with parameter fi = AirnT / B 2 « 1 where 
T is the mean plasma energy, because the jet remains force-free when particles gain small energies in the process of 
acceleration. Parail & Pogutse (1986) show that instabilities may arise from the anomalous Doppler resonance of the 
waves with the accelerated tail of the distribution function. The expression they find for the increment of instability 
is 



7k 



2-k F ™ o J 



dpS^-k^) — fc„ +-Wk|2 



k\p\ 



dpS^ + ^-k^)-^-^ — 



Of 



J B0 



df 



P± dp± 



(16) 



Here lub = Wbo(1 — u 2 /c 2 ) -1 / 2 and / = f{jp\\,Pi_) is the distribution function of particles in drift variables p\\ , p±. The 
first term corresponds to the Cherenkov resonance and the second corresponds to the anomalous Doppler resonance. 
This expression was written with the assumption that k±p±/m c LJBo < 1- Parail & Pogutse (1986) demonstrate the 
instability for the distribution function of electrons accelerated in an electric field and call it the "fan" instability. This 
instability is also present for all distribution functions with small p± and large p^ . We consider a particular distribution 
function, Gaussian in p± and power-law with AN index of a in with the minimal longitudinal momentum p m ; n . 
The resulting 7k is expressed as: 



7T w k w 2 a - 1 ( p 
7k = -rr m, 1 



2 ujbo k 2 p„ 



\uj B om c 



( wbo 
-a 

V ^k 



I 2 , 



1 g 
2fcf 72 



(17) 



where 71, 72 are the 7- factors of the particles in Cherenkov resonance and anomalous Doppler resonance corre- 
spondingly. If 71 w 72 w 1 (the particles in resonance are non-relativistic) we see that instability arises when 



2a 



( UBO 



k 1 m'ui 



2, ,2 

BO 



, or p ± < 



BO 



( 



a+l 



i.e. if p± is small enough. Quasilinear analysis of 



V wk / &| k 2 pj_ ' ' 2a k 2 V / 

the influence of the excited waves on the distribution function shows (Parail & Pogutse, 1986) that the developing 
instability leads to the isotropization of the distribution function. This may easily be seen from the resonance condition 
ftcjk = nhujB + (1 is an integer), which may be treated as the energy conservation law for one interaction of 

a particle with the wave. Here hk\\v\\ « — [(m V|| — ?ifc||) 2 — m 2 v 2 ]/ (2m c ) is the decrement of the longitudinal energy, 
which is positive in the case of anomalous Doppler resonance (n < 0). So the longitudinal energy decreases, and the 
transverse energy increases (n < 0). In the case of Wk << Wb the energy transmitted to the waves is small. Due to this 

CO BO 

instability the distribution function of electrons becomes isotropic with deviations from isotropy of the order of ~ , 

where r is the smallest characteristic time of non-stationary processes increasing the anisotropy of the distribution 
function, and the plasma frequency is w 2 = Aime 2 /m e . Given typical values of B ~ 10 _2 G, n w 0.1 cm -3 we have 
u)bo/u 2 ~ 10 -3 sec. This time is many orders of magnitude smaller than any possible r either due to synchrotron losses 
t s or w or ss l/ft F . Therefore, the distribution of accelerated particles must be highly isotropic in momentum 
space. 

Although the case with > 1 (strong accelerated particles) is apparently contradictory with our assumption of a 
force-free jet, the process of acceleration takes place in the close vicinity of the Alfven resonance and does not affect 
the overall structure of the fields in the whole jet. As we show later the process of acceleration leads to the power-law 
spectrum of particles with high characteristic energy (see Eq.(22)) and the minimum energy a few orders of magnitude 
lower. So the case with (3 > 1 can occur. 

In the case of (3 > 1 there exist many types of instabilities, some of which are hydrodynamic. We will treat plasma 
instabilities according to Mikhaylovskii (1975), who showed that in the case p|| > p±, perturbations with fc|| = have 
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maximal growth rates. The frequencies of these perturbations with k _L B are the solutions of the dispersion relation 
£33 — N 2 = 0, where e 33 is the component of the dielectric tensor along the magnetic field, TV — ck/uo. The dispersion 
relation of the small oscillations has the form 



.21.2 ,.,2 



c 2 k 



Aire 2 Aire 2 



dp < v 



l-J$(kv±/w B )) 



dp\\ u v±_ dp 



00 1 

71=1 B j 

The anisotropic instabilities result in the excitation of potential waves which change the direction of the particle 
momentum p but do not change the energy £ of the particle. Indeed, momentum conservation in the process of 
emitting or absorbing a wave reads 

Sp = hk = Two ( — — 

\ c k 

where N is the wave refractive index. For the potential waves, N » 1. If Sp is not orthogonal to p, then S£ w 
|<5p| • c = huj ■ N » huj. This is in contradiction to the energy conservation S£ — fhv. So, <5p _L p with the accuracy 
1/JV«1. 

We see that the effect of anisotropic instability is to change the particle distribution function in the direction of 
momentum, and not in energy. Because the instability is rapid, the quasi-stationary distribution function is determined 
from the condition 7k = 0. This condition is the equation which should be used instead of Eq. (15) describing the 
change of Jj_ only due to synchrotron losses. The condition 7k — is the integral equation for the particle distribution 
function f(p\\,p±), which follows from Eq. (18). It is difficult to solve this integral equation without knowledge about 
the shape of f(p\\,p±)- We estimate the integral in the dispersion relation (18) assuming the cold distribution function 
/, although in reality it is smooth in momentum space. Also, we assume that kv± is of the order of lob and estimate 
expression (18) taking into account only the first terms of The expansion in Bcsscl functions over its arguments. Then, 
the threshold of stability (7k = 0) when particles have equal p± and p» in the ultrarelativistic case is found to be 

P\\ _ 4-/3+^(17/3-8) 

p\- 4(/?-l) ' liyj 

where = Airnm c c 2 -f / B 2 . We denote the threshold of stability p\\/p±_ = a, a > 1. Expression (19) is not quite 
valid for very large /3(> 8) but large (3 are not achieved here. 7k > arises when pu > ap±. The instability is 
fast because its growth rate is proportional to the cyclotron frequency; still the growth rate is much less than the 
cyclotron frequency. We will assume that the process of acceleration and the fast instability considered above will form 
a nonequilibrium distribution function with anisotropy of order unity. However, after particles leave the acceleration 
region (/3 > 1), various instabilities such as the "fan" instability considered in this section lead to fast full isotropization 
of the distribution function. 

Now we understand that dJ± /dt ^ not only because of the synchrotron losses but also because of the anisotropic 
plasma instability. We use the relation p\\/p± ~ a as an approximate closure equation for the system of two equa- 
tions (1), and express the longitudinal and perpendicular components of the particle momentum through the energy 
by the relations 

p2± = (1 + a 2 ]c 2 5 P " = (1 + a 2 )c 2 

Since the dependence of a on (3 (equation (19)) is weak for > 2, the value of a is approximated as a numerical 
constant independent of r in further expressions for particles trajectories and particles distribution function. For the 
derivative of the energy we write equation (11) adding the synchrotron losses. 

d£ _ £v\ (lo + k)[F - uj] dx 2p\ 



£ 2 - mlc 4 2 a 2 (£ 2 — mlc 4 ) 



dt 1 - v|, rS dt 3 TO C 50 c 

Substituting dx/dt from (14) we obtain 

d£ _ £ {uj + k)[F- lu] /wMiBi\ 1/2 2 £ 



u 2 m r e . (20) 



,2 



dt l + ^-aVTTa 2 rS \ At J 31+a 2WB0 m e 

Let us introduce some notations for convenience: 

1 (u + k)[F - u] 



(21) 



1 + a 2 — a\J\ + a 2 
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3 m c f(w + k)[F-w}\ uj* 2 A x B 



1 4 1 + 2a 2 - 2aVTT^ \ S J u 2 B0 r c r 2 uj* ' 

Here (ui + k)[F — w]/S is fixed for given wavenumbers of the perturbation, m and k. It is large near the fast magnc- 
tosonic resonance surface 5 = 0. For an AGN jet where B 10~ 2 G (Begelman at al., 1984), r w 1 pc, the coefficient 
l/(ujg r c r 2 LU*) i=a 10 4 , which is a large number. The value of uj* 2 A\B\ is the dimensionless amplitude of the per- 
turbation squared as we noted after Eq. (13). Now taking into account our notations, the solution for £ has the 
form 



(23) 



Here £o = £{x = 0). Let us denote Qx/r = x' . If x' is of order unity, then £ ~ £q and particle acceleration is 
not effective. We assume that x' can be large because the quantity Q is large in the case of maximal damping rate 
of eigenmodes. If x' is large, then £ ^> £ q and the particle acceleration is effective. It is seen that first £ increases 
exponentially and then decreases as 1/x'. We consider the asymptotic behaviour of £ in the case of different x' and 
£ o . In the case of £ < £\ 

£i/x', x'e x >£i/£ ^ 
£ e x ', x'e x ' <£i/£ ' 



In the case of £q > £\ 

£ = 



£i/x', x' > 1 



2^/x' 2 , x' < 1 



We ignore the case when £q > £\ since the typical initial energy is much less then £\. Knowing the particle trajectories 
(Eqs. (14) and (23)), we can find the distribution function F(£,x) of accelerated particles in the phase space of (£,x) 
using the fact that the trajectories of the particles are the integrals of the collisionlcss kinetic equation. 



3. Formation of the spectrum of accelerated particles 

Our goal is to calculate the distribution function averaged over x knowing the distribution function at x = 0. We 
will consider the stationary case, when df/dt = 0. Let the distribution function be given at point x' as a function of 
£' = £ /£ i . Then the number of particles 

dN = F Q (£ Q ,x' )d£ Q dx' = F (£' Q {£',x'),x' {£',x>)) ^^j d£'dx'. 

Below we scale x with r/Q and £ with £\ and omit primes in the dimensionless x' and £' . Thus if wc know the trajectory 
£o{£, x),Xq(£, x) we may obtain the distribution function at the point £, x. In our case the Jacobian D(£q, xq)/D(£, x) 
is not equal to unity because the variables £,x are not canonical. So, 

F(£, x) = F (£ (£, x),x (£, x)) ^°'^ ) . (25) 

If we repeat the derivation of Eq. (23) with an arbitrary reference point x n ^ (£(x ) = £ ) we obtain 

£ = {x-l+ (^^y -xo + lj cxp(-x + xo)}- 1 . (26) 

Substituting £ expressed from (26) into (25) we obtain 

m „ - « ({-i +I0 + (I - * + 1) ox P „ - (§^ - . m 

F ^|-1 + + - x + 1^ cxp(x - xo) | ,X(^j ^| cxp(x — xq)^-. 

We know that initial fast particles are injected near the point x = . Therefore in the last formula we let xq = 
everywhere except in the denominator. Physically F (£o, 0) is the particle density at zero point and x is the distance 
below which our formulae, for instance (14), are not valid. Because of the small imaginary part of w which is an 
attenuation increment of the wave (Istomin & Pariev, 1996), the divergence In a; at the resonance is cut at the point 
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where Imw <C Rew. Here A(r — r\) is not equal to zero as we assumed earlier, but determined by the small value of 
the imaginary part of u. Thus, x is the width of the resonance x = Q(ro — r A ) / r A ■ Eventually, Eq. (27) is 



F(£,x)=F - x + l) exp(x)} ,x ^j 



||exp(x)-. (29) 
£ z x 



Because the width of the resonance Xq and the width of the acceleration region Ar , which is some function of xq , are 
small compared with the jet radius we average the distribution function over the jet radius near the resonance surface 

oo 

F(E) = / F(£,x)dx, Ax = —Ar. (30) 

Ax J TA 

o 

Denoting N(£) = F(£)Ax we obtain 

N(£) = J F U-l+(^-x + ljexp(x)} ,x ) ^exp(x)^-dx. (31) 
o \ / 

Let us change the variable of integration to £ - The dependence of £q versus x under fixed £ is presented in Fig. 1. 
The value of £ m i n is 

£ mi „ = l/(-l + exp(l/£)). (32) 
Because the function x = x(£o) has two branches, the equation (31) is transformed into 

£ 



^)~/ F (£ ,x Q)Y ^- xd £ 0+ 



10 



A.R. Beresnyak et al.: Relativistic parsec-scale jets: I 



oo 

/ Fq(£ ,xo)—— d£ , (33) 

J It - x 



t Z X J I/O 

£mi„ 

where x is taken on the first branch of x(£q) in the first integral and on the second branch in the second integral. Let 
us make some natural assumptions. The initial distribution function is cut off at some energy of T , provided that 
T <C 1 (bearing in mind that T as well as £ are in the units of energy £\ (22)). Let us calculate N(£) in the region, 
where T <C £ <C 1, because it is seen from Eqs. (32) and (33) that in the case £ ~ 1 the integral is small. Let us 
evaluate the contribution of the singularity x = l/£ in integrals (33) 

p _ p _ d 2 £ a ( 1 

£-0 £rnin — 7 9 I *^ r> 

The integral near the singularity is 



dx 2 V £ 



1 . \\jirr 

"?2 / ^b(£minj x ) - d£ ~ \/£o — £ n 



Thus, its contribution is small. Because the main contribution to Eq. (26) is when x is far from l/£, we may use 
asymptotic values of x (see Eq. (24)): x rts ln(£/£ ) on the first branch; x w l/£ + 1 on the second branch and obtain 

£ oc 

fmin £mi„ 

Let us consider only £ min <C T i.e. £ <C — 1/ In To because in the reverse case the integral is small. Now we can replace 
£min by and also l/£ 3> ln(£/£ )- Eq. (34) becomes 

£ 

N(£) = -^- J F a (£o)(-\n£o)d£ + (35) 
o 

where n = J °° Fq(£q) d£$. Rewriting Eq. (35) as averaging < . . . > over the initial distribution function Fq 

"< f » = 5 

we note that the first term is approximately equal to 



N ^ = £x- <-^> + ^ ^ 



7^ (-In To), 
£x 

and is small compared to the second term because the inequality l/£ 2 ^> — lnT is the consequence of inequality 
l/£ ^> — lnT . Therefore N(£) cx l/£ 3 , and we obtain a power-law distribution with the index —3. The validity region 
of this formula T < £ < — 1/ lnT is quite large. 

We performed the integration in Eq. (33) numerically, with Gaussian initial distributions at different temperatures. 
The results are presented in Figs. 2 and 3. 

The value of the power-law index is —3 for small energies and rises slightly before the cutoff of the distribution 
function. However, the —3 power-law index is not applicable all the way down to £ = 0. The natural boundary is the 
width of the acceleration region, Ar (see Eq. (33)). Because the main contribution to the integral (28) is from the 
region near the point x — l/£ + 1, the lower limit for power law l/£ 3 is £\ = 1/Ax. We also found the mean power 
law index for the distribution of particles between energies £\ and — 1/lnTo. The results are presented in Fig. 4. 

Recall that in the acceleration region the distribution function of the accelerated particles has an anisotropy of 
order unity, but the distribution of particles escaping the acceleration region becomes isotropic in momentum space 
very quickly. 

One remaining question in the formation of the spectrum of energetic particles is the origin of the initial particles, 
since to first order in the perturbed fields there is no acceleration of the cold plasma in Eqs. (1). We have already 
mentioned that the particles to be accelerated originate in the region close to the Alfven resonance. It can be seen 
from Eqs. (1) that to second order in the perturbed fields, acceleration always takes place and the power acquired 
corresponds to the term j • E in the energy conservation equation for the electromagnetic field. The damping of the 
wave occurs due to the heating of the plasma by current ji which is described by second-order terms. If all of the wave 
energy is transmitted to the particles, the condition f3 > 1 would imply that the energy density of the wave is greater 
than the energy density of the stationary field. Let us assume this at least on the surface of the resonance where the 
field is maximal. Then the energy acquired by one particle can be evaluated as T w Bq/Sttu. Given typical values 
B a w 1CT 2 G and n w 0.1 cur 3 we obtain T w lOMeV. 
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Fig. 2. Distribution function calculated with a Gaussian initial distribution with temperature 2 • 10 5 £i (the energy 
£ is in units of £\ (22)) 
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Fig. 3. Distribution function calculated with Gaussian initial distribution with temperature 10 3 (the energy £ is in 
units of £ x (22)) 
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Fig. 4. Averaged power law index (the initial particle energy T is in units of £\ (22)) 



4. Summary 

In the present work we considered possible acceleration of electrons and positrons inside relativistic rotating electron- 
positron force-free cylindrical jets with spiral magnetic fields. It is very plausible that inner parsec-scale jets in active 
galaxies can be described in the frame of this physical model. The observations of jets in 3C273, M87, 3C279, 3C84 
and PKS0528+134 show that extragalactic jets are likely to be electron-positron ( Morrison & Sadun, 1992; Reynolds 
et al., 1996; Wardle et al., 1998 ). It also seems that such jets are not in steady state and can contain waves excited 
by environmental processes, of which the most powerful would be the variability of the accretion rate onto the central 
black hole and accretion disc instabilities (magnetorotational instability). Both strongly perturb the magnetic field 
at the base of the jet. We considered the behaviour of such excitations inside cylindrical force-free jets embedded 
in a dense medium in our previous studies (Istomin & Pariev, 1994; Istomin & Pariev, 1996). We estimate that the 
density of the gas commonly found in the nuclei of galaxies is large enough that it should be taken into account when 
considering the evolution of perturbations inside the force-free magnetized jet. We found that in the wide range of 
parameters determining the equilibrium structure of electromagnetic fields and perturbations inside the jet there exist 
resonant surfaces on which the phase speed of eigenmodes is equal to the local Alfven velocity. Those eigenmodes for 
which there is an Alfven resonance have a small damping rate, whilst those having no Alfven resonances are neutrally 
stable. This shows that the existence of an Alfven resonance leads to losses of energy of corresponding eigenmodes. 

Indeed, the amplitude of the wave is large in the vicinity of the Alfven resonant surface and decreases away from the 
Alfven resonant surface. Both electrons and positrons are subject to drifting out of this region of strong electric field of 
the wave due to the non-uniformity of the electric field. When moving to a region of weak electric field, particles gain 
energy. This energy is taken from the energy of the wave, which is an eigenmode in the jet, and the eigenmodc decays. 
In this work we propose a novel mechanism for acceleration of particles in situ in strong electromagnetic waves of 
varying amplitude inside the jet. Unlike the acceleration on shock fronts proposed to operate inside bright knots (e.g., 
Blandford & Eichlcr 1987), the acceleration by electromagnetic waves is regular, i.e. it does not require turbulence or 
inhomogeneities to be present in the jet. Also, the net gain of energy by particles in our mechanism occurs slowly on 
a time scale greatly exceeding the Larmour orbital time and also exceeding the period of the wave. In this aspect, the 
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regular electromagnetic acceleration by a nonuniform wave differs from acceleration in magnetic reconnection layers 
(Romanova & Lovelace 1992; Blackman 1996), where strong departures from ideal MHD occur and a particle gains its 
total energy during one passage through the region of strong electric field. Our acceleration mechanism does not rely 
upon violation of E • B = in any region inside the jet as the gain of energy during the drift is a result of a kinetic 
treatment of plasma possible only when one takes into account a finite Larmour radius. 

We found how particles are accelerated in the strong fields near the Alfven surface. We considered the process 
of particle acceleration using equations of motion in drift approximation with addition of synchrotron losses and 
isotropization of particle distribution by plasma instabilities. Such instabilities are excited by accelerated particles 
if the distribution function is anisotropic in momentum space, and act to conceal any anisotropy of the distribution 
of particles. Acceleration process and synchrotron losses taken together form a power law energy spectrum of ultra- 
relativistic electrons and positrons with index between 2 and 3 depending upon the initial energy of the injected 
particles. This is consistent with the typical spectral indices of radio emission observed in parsec-scale jets. The 
power law spectrum extends up to the energy £ max , where a sharp cutoff occurs (see Figs. 2,3). The magnitude 
of £ max depends on the initial particle temperature To as well as on the characteristic acceleration energy E\ (26): 
£max ~ £i/ln(£i/T ). The quantity T evaluated from the equipartition condition in the acceleration region is equal 
to 10 McV for the magnetic field 1CP 2 G and density 0.1 cm~ 3 correspondingly. The quantity £\ is, according to (26), 
of the order of 10 4 (5B/B ) 2 MeV for the value of B a = 10~ 2 G, where 5B/B is the dimensionless amplitude of the 
perturbation. For large perturbations (SB w B ) £\ ~ 10 4 McV. Thus particles accelerated near Alfven resonance are 
in the energy range of 10 MeV < £ < 10 3 MeV. These particles emit synchrotron radiation in the range of frequencies 
from approximately 100 MHz to 1000 GHz which covers the frequencies of modern radio observations. We present 
results of the calculations of synchrotron emission and polarization produced by the accelerated particles with power 
law distribution over energy in a separate article (Pariev et al. 2002). 
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